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We study the self force acting on a scalar charge in uniform circular motion around a Schwarzschild 
black hole. The analysis is based on a direct calculation of the self force via mode decomposition, 
and on a regularization procedure based on Ori's mode-sum regularization prescription. We find 
the four self-force at arbitrary radii and angular velocities (both geodesic and non-geodesic), in 
particular near the black hole, where general-relativistic effects are strongest, and for fast motion. 
We find the radial component of the self force to be repulsive or attractive, depending on the orbit. 
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PACS number(s): 04.25-g, 04.70.-s, 04.70.Bw 

The problem of finding the equations of motion for 
a particle in curved spacetime has become recently ex- 
tremely important, as the first generation of interfero- 
metric gravitational wave detectors will soon be opera- 
tional, and with the prospects of having a gravitational 
wave space antenna in the not-very-distant future. The 
generation of very accurate templates for the waveforms 
detected from a system of a compact object in orbit 
around a supermassive black hole is an extremely hard 
task. It is likely that one would need to have accurate 
templates for as many as 5 x 10 5 orbits. For such a 
system, accurate templates are necessary for detection, 
because the predicted signal-to-noise ratio for LISA is 
approximately of order 10 for a 1 year integration time. 
Lack of accurate templates would result in a loss of a fac- 
tor of roughly the square root of the number of orbits in 
sensitivity jjj , which would result in signal-to-noise ratio 
below the detectibility threshold. 

In order to generate accurate templates, an important 
necessary ingredient is the inclusion of radiation reac- 
tion in the orbital evolution of the compact object. The 
radiation-reaction forces need to be calculated locally, 
i.e., in the neighborhood of the compact object. In the 
conventional approach, one calculates, at infinity and at 
the event horizon of the black hole, the fluxes of quanti- 
ties, which are constants of motion in the absence of ra- 
diation reaction. Then, one uses a balance argument to 
relate these fluxes to the rate of change of a correspond- 
ing local quantity of the compact object. This approach 
generally fails because of the inadditivity of the Carter 
constant. For very simple cases, e.g. for circular or equa- 
torial orbits around a Kerr black hole, the evolution of 
the Carter constant is trivial, such that the conventional 
approach is useful [||. However, generic orbits around 
a rotating black hole arc neither circular nor equatorial, 
and consequently a new approach, which is not based on 
balance arguments, is of great need. 

Several approaches have been suggested for the calcula- 
tion of self forces. Quinn and Wald || and Mino, Sasaki, 
and Tanaka j| recently proposed general approaches for 
the calculation of self forces. However, it is not presently 
clear how to apply these approaches directly for actual 



computations, the greatest problem being the calculation 
of the so-called "tail" term of the Green's function, which 
arises from the failure of the Huygens principle in curved 
spacetime. (In addition, there are in general also lo- 
cal, Ricci-curvature coupled and Abraham-Lorentz-Dirac 
(ALD) type terms ||, which are much easier to calculate, 
the former vanishing identically in vacuum.) 

Recently, Ori proposed a local, causal approach for 
the calculation of the self forces which is based 

on decomposition of the self force into modes, and on a 
mode-sum regularization prescription (MSRP). Although 
MSRP is not fully proven as yet, it has already been 
shown to be valid for simple cases, such as scalar charges 
in general orbits in Schwarzschild spacetime, and, in par- 
ticular, for circular orbits which we consider here. MSRP 
is likely to be susceptible of generalization also for mas- 
sive particles in orbit around a Kerr black hole. If ro- 
bust, MSRP can be of great importance for the gener- 
ation of accurate templates. We hope that MSRP can 
be combined with other approaches, which were recently 
proposed, such as mode decomposition of the self forces 
which are sourced by just the distant past world line 
or a normal-neighborhood expansion JTo| . 

We first describe very briefly the main ideas of MSRP 
{7|,|| . Then, we apply MSRP for the case of a scalar par- 
ticle in circular orbit around a Schwarzschild black hole, 
and calculate the self four-force acting on the particle 
linearized in its own self field. 

The contribution to the physical self force from the 
tail part of the Green's function can be decomposed into 
stationary Teukolsky modes, and then summed over the 
frequencies uj and the azimuthal numbers m. The self 
force equals then the limit e — > CP of the sum over all 
I modes, of the difference between the force sourced by 
the entire world line (the bare force bare i^) and the force 
sourced by the half-infinite world line to the future of e, 
where the particle has proper time r = 0, and e is an 
event along the past (r < 0) world line. Next, we seek a 
regularization function /if, which is independent of e, such 

that the series J2i( haie Fft ~ ^) converges. Once such a 
function is found, the regularized self force is then given 
by ™ n F, = £,( barc i^ - hi) 



where is a finite 
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valued function. MSRP [|7|,|| then shows, from a local in- 
tegration of the Green's function, that the regularization 
function = a^i+b^+c^i -1 , and for the case of a scalar 
charge in circular orbit around a Schwarzschild black hole 
MSRP yields the values of the functions a M , b^, c M and 
analytically. In particular, it can be shown that for such 
orbits = = c M and d r — 0, such that in practice the 
regularization prescription of the radial force is reduced 
to subtracting b r from each £ mode of the bare radial 
force. Note that & M is just the limit I -> oo of barc i^. 

In the following we describe the results obtained from 
this new approach for the case of a point-like scalar 
charge in circular orbit around a Schwarzschild black 
hole. Our results are fully relativistic, i.e., we do not 
introduce any simplifying assumptions such as far field 
or slow motion. Because of the fully relativistic na- 
ture of this study, our analysis is numerical. How- 
ever, it is reasonable to expect that analytical solutions 
will not be available in general, except, possibly, only 
for very simple cases, such as static configurations [ ]lT| . 
We calculate the contribution to the force which the 
scalar charge feels, due to its own field, to leading or- 
der in the particle's charge. We use spherical Regge- 
Wheeler coordinates, for which the Schwarzschild metric 
is ds 2 = (1 - 2M) (- di 2 + dr* 2 ) +r 2 ( d0 2 + sin 2 9 d<j> 2 ) , 
where M is the black hole's mass, and the radial 
Schwarzschild coordinate r(r*) is given implicitly by 
r* = r + 2M\n\r/{2M) - 1|. The field satisfies the 
wave equation V li 'V ti ^(x a ) = — 4irp(x a ), where the 
charge density p(x a ) = q drS 4 [x a — Xq (t)](— g)^ 1 ^ 2 - 
Here, q is the charge of the particle, whose world line 
is Xq(t), t being its proper time, g being the met- 
ric determinant, and V M denotes covariant differentia- 
tion. We take the charge to be in circular orbit at 
r* = rj), = ir/2, and d(j)/ dt = il. (We are not 
restricted to the Keplerian angular velocity £Ik-) We 
next decompose the scalar field $ into modes according 
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where W& m is the Wronskian determinant of ^£ m (r*) and 
*L(r% T em (r*) = ^} m (r*m m (r*), and S lm (r*) = 

*tm(r*)*7m,A r *) + nm(r*)*L,A r *)- We find that 
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We also obtain baro i^™ = _0- lbare F/ 



and ^re pirn = 

0. It is convenient to define new radial functions 

Z ±£m( r *j by ^±to( r *) = giwr'^ite^^ which gat _ 

isfy the homogeneous equations 



d 2 Z ±lm dZ ±lm r S1 
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with boundary conditions Z + lm {r* > M) = 1 + aff+ + 
4fl + O(fl) and Z- lm {r* < -M) = 1 + a^/- + 
a^fl + O(fl), where /+ = (wr) _l , /_ = 1 - 2M/r, and 
a+ = i£(£+l)/2, a+ = [-£(£- l)(£+l)(£+2) + 4icoM}/8, 
<q = [1 + £(£ + 1)]/(1 - 4iwM), and = [£(£ + 
l){£ 2 + I + 6) + UujM + 4]/[4(l - 2iujM)(l - tewM)}. 
We solve Eqs. (^) numerically using both Burlisch-Stoer 
and fourth-order Runge-Kutta integrations with adap- 
tive step-size controls. Both integrators yield results 
compatible within the error limits. We place the ex- 
terior and interior boundaries at a distance of several 
to $ = du>J2e m e ~ l,Jjt ^£™(r*)Y ern (9,(f))/r(r*), such mode wavelengths (in r*) from rj, and then use sue- 
that the equation which governs the field ^ em becomes cessive Richardson extrapolations, with increasing dis- 
tance to the boundaries, until the extrapolation of the 
boundaries to r* — > ±oo yields an error smaller than 
a given threshold. Notice that for modes with m = 
the wavelength is infinite, such that the boundaries can- 
not be taken far enough from the charge. Indeed, we 
find that for this case the Richardson extrapolations do 
not converge. Instead, we can solve for this case an- 
alytically. We find that f^ m=0 = and f^ m=0 = 
(2TT/ 1 )(q/M) 2 Y eo2 (7r/2,0)[2Q e (p )dP e (p a )/dp+ 1/(1 - 
Pq)], where p = (r — M)/M . Here, Pi, Qi are the Legen- 
dre functions of the first and second kinds, respectively. 
Figure [l] shows the functions Z^ =1 m=1 (r*). Similar qual- 
itative behavior is found also for the other modes. Until 
the peak of the effective potential barrier the functions 
Z± im var y on iy slowly, and then start oscillating rapidly. 

The temporal component of the bare force is finite. 
(In the regularization scheme this corresponds to b t — 
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This equation should be solved for each mode £m 
with boundary conditions of ingoing waves at the event 
horizon (r* — > — oo), and outgoing waves at infin- 
ity (r* — > oo). The effective potential is given by 
Vi(r ) = (1 - 2M/r) [2M/r 3 + £(£ + l)/r 2 ], and 7 = 
1 / — 2M/i — r 2 2 . The contribution of the £m mode 



to the force is given by barc i^™ = g(<j>< 



(note that for circular orbits u"& e ™ = 0). 

We next find numerically the solutions ^tmi^Jm) f° r 
the homogeneous equations corresponding to Eq. (Q), 
which satisfy the boundary condition at infinity (the hori- 
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0.) MSRP predicts that dt exactly balances the ALD 
force, such that the full self force is given only by F t = 
bare F t . We compare our results with their Minkowski 
spacetime counterparts. In flat spacetime one can solve 
analytically for each mode, sum over all modes, and find 
that F t Min = Iq 2 ^ 2 ^^. Here, 7 Min is the usual flat 
spacetime Lorentz factor. For the comparison we choose 
the same values of ro , for the curved and flat space- 
times. Figure || displays Ft as a function of r/M for two 
cases: (|]A) Non-geodesic circular orbits, with a fixed an- 
gular velocity Ql. (In this case the tangential velocity 
increases linearly with r.) When r/M is large the value 
of F t approaches its flat spacetime counterpart. Recall 
that il = d<f>/ dt, where t is the time of an observer at 
infinity. Because of the red-shift effect at small values of 
r, orbits with the same value of SI have very large proper 
tangential velocities at small r/M: a fixed SI corresponds 
to the ultra-relativistic limit when the orbit is close to 
the black hole. The second case is (^B) geodesic motion, 
which satisfies Kepler's law Sl^r 3 = M. The innermost 
(unstable) causal orbit is located at r = 3M. Approach- 
ing r = 3M, the motion of the particle approaches the 
ultra-relativistic limit, which is manifested by the rapid 
growth of Ft. At larger radii the value of F t approaches 
the flat spacetime counterpart. This is detailed in Fig. 
|(C). 





FIG. 1. The radial functions Z ± for I = 1 and m = 1 as 
functions of r*/M. The particle is at = 4M and = Qk- 
Top panel: Re[Z + (r*)] (solid line) and Re[Z~(r,)] (dotted 
line). Bottom panel: Same as in the top panel, for Im[Z ± (r»)] 

Next, we study the radial, conservative component of 
the self force. First, we check the agreement of our 
numerical results with MSRP. MSRP predicts f| that 
a r = = c r and that 
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We check the accuracy of our numerical determination 
of the value of b r by comparison to Eq. (if) . This check 
serves the two-fold purpose of (i) checking the numerical 
code, and (ii) checking the compatibility of the analyt- 
ical prediction of MSRP for the regularization function 
h r with the numerical determination of the bare force. 
Figure |(A) shows the behavior of barc i^, after summa- 
tion over m and w, as a function of £, and Fig. ||(B) 
displays barc Fjf — 6 r as a function of I. This difference 
decreases like t~ 2 for large values of £, which confirms 
the predictions of MSRP for the values of a r , b ri and c r . 
(We emphasize that we cannot test numerically the pre- 
diction of MSRP that 4 = 0.) Note that the radial ALD 
force vanishes, such that Yen F r is the full self force. Sim- 
ilar behavior is found also for Keplerian orbits at other 
values of 7q, and also for non-geodesic circular motion, 
with angular velocities both greater or smaller than Qk- 




FIG. 2. The temporal component of the self force Ft as 
a function of r/M. Top panel (A): Fixed orbital valocity 
Q. — 3.2 x 10 _2 M _1 (non-geodesic orbits). Q = Qk when 
ro = 9.92M. Bottom panel (B): Keplerian orbits. Circu- 
lar orbits in Schwarzschild are depicted by solid line, and in 
Minkowski by a dashed line. Bottom panel (C): The rela- 
tive difference of the Schwarzschild and Minkowski results for 
Keplerian motion. 

The regularized component of the radial self-force 
tall F/ is obtained by subtracting b r = °axepe->oo f rom 
each mode harc F^. The total regularized force is then 
obtained by 



prcn \ ^ tail 

r ~~ i r 
n=0 



■11 



e+i 



(6) 



(5) 



where the remainder lZf, +1 is given approximately by 
n i+i _ £2(bare F) « _ 5 r )^(i)(£ + i) f or sufficiently large 

values of £, ^ W (x) = d 2 In T(x) / dx 2 being the trigamma 
function. Note that for large arguments %p^(x) w x^ 1 . 
Figure p| displays the regularized radial self-force for Ke- 
plerian (|4]A) and non-Keplerian (|^B) orbits. The radial 
self force in the far field limit (r 3> M) is repulsive, and 
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satisfies F™ n w aq 2 M 3 r~ 5 , where a is a dimensionless 
parameter of order unity. A minimum-^ 2 fit shows the 
exponent of r to equal —5 and found a to equal unity, 
both with 3% errors. However, in the strong field, the 
force law deviates from this simple relation, and grows 
faster. For the non-Keplerian orbits, we find that in the 
slow motion limit (f2 <C Qk) the radial force is propor- 
tional to f2 2 . The exponent of fl is found to be 2 with a 
3% error. Combined with the result for Keplerian orbits, 
we find that for any circular orbit, in the far field and 
slow motion limits, the radial force is repulsive, and is 
given by 



f 1 ': 



aq 2 (G 3 /c 6 )M 2 n 2 /r 2 
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This results explains the vanishing self force in the static 
limit fljjllj . For faster motion the fi 2 law does not hold 
any more. In fact, for £1 > Hk we find that the radial 
self force varies rapidly, and eventually changes from re- 
pulsive to attractive. The radial self force does not cause 
a net change in the energy of the particle. However, if 
the orbit has a non-zero eccentricity, this force induces an 
additional precession of the periastron, which in the slow 
motion or far field limits is retrograde. This precession 
has an effect on the frequencies of the emitted radiation. 



r/M 



FIG. 4. The regularized radial self force. Top panel (A): 
F r rcn as a function of r/M for Keplerian (geodesic) orbits. 
Bottom panel (B): F r Tcn as a function of fi/fix for rg = 4M. 
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FIG. 3. The radial component of the I multipole of the 
bare force as a function of I. For these data the orbit is at 
ro = 3M. (A): baro F* for a geodesic orbit. (B) barc F r £ - b r for 
the same orbit. 

Our results show that the self force can be calculated 
for a simple, although non-trivial, problem. MSRP was 
found to be useful also for other cases, e.g., static scalar 
and electric charges in Schwarzschild for scalar and 
electric charges in circular motion in flat spacetime [p"3| , 
and for general radial motion of scalar charges in spher- 
ical symmetry jl4j . A closely related approach was used 
also for a mass point in radial free fall in Schwarzschild 
]15| . We hope that similar methods can be used for more 
realistic cases, which may be relevant for the orbital evo- 
lution of compact objects around black holes. 
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